Pattern formation during diffusion limited transformations in solids 
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We develop a description of diffusion limited growth in solid-solid transformations, which are 
strongly influenced by elastic effects. Density differences and structural transformations provoke 
QQ , stresses at interfaces, which affect the phase equilibrium conditions. We formulate equations for 

■ the interface kinetics similar to dendritic growth and study the growth of a stable phase from a 

' metastable solid in both a channel geometry and in free space. We perform sharp interface calcu- 

, lations based on Green's function methods and phase field simulations, supplemented by analytical 

investigations. For pure dilatational transformations we find a single growing finger with symmetry 
breaking at higher driving forces, whereas for shear transformations the emergence of twin struc- 
) ^ tures can be favorable. We predict the steady state shapes and propagation velocities, which can 

^— ( . be higher than in conventional dendritic growth. 

I. INTRODUCTION 

^ ■ Transformations between different solid states of a material are essential for many technological and scientific appli- 
cations, and understanding their kinetics is therefore not only interesting because of the scientific variety and beauty, 
but also important for tailoring of new materials with specific properties. Often, these processes are accompanied 
by elastic deformations, which can be e.g. due to density differences or structural changes, which provoke stresses 
at interfaces between adjacent phases. The influence on the thermodynamics of the transitions between different 
phases has been thoroughly discussed in the literature [l|, Q , whereas the kinetics of these processes are still far less 
understood. 

Facing this lack of understanding, we would like to mention that very fruitful concepts have been developed to 
describe the dynamics of solidification or melting processes as moving boundary problems. It means that the evolution 
^ of the different phases is described on the very detailed level of the motion of the interfaces between them, for which 
Q ' suitable equations of motion have to be posed. A particular application is the famous example of dendritic growth, 
O . where heat or component diffusion is the rate limiting process. At the interface between the solid and melt phase the 
temperature or concentration is close to its phase equilibrium value. The concentration jump or the release of latent 
, heat due to partitioning leads to a jump in the corresponding fluxes at the interfaces, and the magnitude of this jump 
J> • is directly related to the local interface velocity. Altogether, this provides a closed description for the motion of each 
CO [ interface point, and many efforts have been done towards analytical and numerical solutions of this problem. 

• It turns out that the questions concerning existence, shape and growth velocity of steady state patterns in solid- 
ification crucially depend on selection mechanisms. In this spirit the influence of many different physical effects on 
^; . solidification or melting processes has been studied throughout the years [1, 0]- For example, the effect of isotropic 
' surface tension was proven not to serve as a selection mechanism for a solution, whereas anisotropic surface tension 
leads to a unique solution. Also, the appearance of triple junctions can provide a selection mechanism, and can 
therefore lead to a solution of the free boundary problem [5i] . 

Of course, the situation considering solid-solid transformation kinetics differs significantly from the solidification or 
melting problem, where elastic effects often play only a minor role. Still, we follow here the same successful concepts 
to formulate a free boundary description for the transformation kinetics in solids, now including the important effect 
of elasticity. In real systems many other effects apart from elastic effects come into play - like inhomogeneous 
cd compositions, the MuUins-Sekerka instability, crystal anisotropics, polycrystalline structures, etc. - and also quite 
successful attempts have been made to model this complex behavior (see for example in Q and references therein). 
Here, in contrast, we focus on gaining a fundamental understanding of the elastic influence on pattern formation 
processes, and we restrict our investigations therefore to simple and "clean" model systems. Correspondingly, we 
consider only a single diffusion field, namely the thermal diffusion in pure materials. 

The formulation of such a model system requires to write down explicit expressions for the boundary conditions 
and suitable equations of motion at the propagating interfaces, which therefore have to be tracked during the entire 
dynamical process. Finally, it provides a detailed description of the microstructure evolution on the level of individual 
grains in a polycrystalline structure. As for the dendritic growth, we focus here on the regime of diffusion limited 
growth. Then, propagation of a front generates or consumes latent heat, which has to diffuse through the system, 
and we assume that this is the rate-limiting process. This is in contrast to isothermal solid-solid transformations that 
were discussed in [X], which are driven by local interface kinetics, therefore allowing fast front propagation also on 
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the scale of the sound speed. Many important properties result from the influence of elastic effects due to structural 
transformations and density differences. 

Apart from sharp interface methods, that are based on Green's function techniques, we also use phase field de- 
scriptions, which are nowadays widely used. More importantly, however, is that we supplement these numerical 
investigations by analytical calculations which allow to understand these processes and the relevant parameters on a 
deeper basis. 



II. FORMULATION OF THE PROBLEM 



Specifically, we consider the growth of an equilibrium solid phase (3 from another metastable solid phase a, as 
depicted in Fig. [1] At a temperature Tgg, both stress-free phases are in bulk equilibrium with each other, whereas 
below this value the new (3 phase is thermodynamically favorable. Far ahead of the propagating front, the system is 
therefore exposed to a lower temperature Too, and we define a dimensionless temperature w = C{T — Toc,)/L, where 
L denotes the latent heat of the transition and C the heat capacity. Within the solid phases the evolution of the 
temperature field obeys a diffusion equation, 

DV^^ (1) 

where D is the thermal diffusivity, which we assume to be equal in both phases, as well as the heat capacity (sym- 
metrical model). 

The motion of the interfaces follows from energy conservation principles, which account here for the net heat flux 
and the generation or absorption of latent heat due to the local motion of the interface with normal velocity Vn , 

Vn^Dn-{S/w^\-Vwt\), (2) 

where n is the interface normal and the superscripts denote the phases. 

Elastic effects enter into the problem via a modification of the local phase equilibrium condition at the moving 
interfaces. In contrast to previous investigations |1, Q, we now assume the phases to be coherently connected, i.e. we 
neglect the appearance of any slips, detachments or defect formation at the interfaces, which is reasonable for small 
misfits. It is important to mention here, that with the assumption of coherency we formulate a completely different 
model being capable to describe qualitatively different types of elastic effects. For example the well-known effect of 
elastic hysteresis p], , which was proven to disappear without coherency Q , is now included in the model. 

In a Lagrangian description of linear elasticity, the displacement field is denoted by u, and the coherency condition 
reads, u^"! — ut'^l at the interface. The strain field is then given by the symmetrical spatial derivative, e^^ — 
{dui/dxk + duk/dxi) /2. This implies that tangential and shear strains eTT,^sT,f-ss are continuous at the interface 
with T, s being the two tangential directions. Force balance demands the continuity of the normal and shear stresses, 
o'nm o'nr, ^ns , where 71 Is the normal direction. Here, the stresses are defined as the derivative of the free energy with 
respect to the strains, (Tik = dF/dcik- 

The interface temperature for a planar, stress free interface is given by = A, with A — C(Teq ^Too)/L. Capillary 
corrections induce a curvature dependent term, and for simplicity we discuss only the case of isotropic surface tension 
7. Finally, elastic deformations lead to an additional change of the interface temperature. Altogether, this leads to 
the local equilibrium boundary condition 10] 

T C ~ 

W^nt = A- doK+ -j^SFel, (3) 

with do = ^Tf.qC/LP' being the capillary length and k the local curvature of the interface, counted positive for a 
convex phase {3 (see Fig. [T]); Notice that due to the coherency constraint [HI, [13], the elastic shift of the equilibrium 
temperature is proportional to the difference of a new potential 5Ff,i = — Fjf' , which is related to the elastic free 
energy i^e/ via the Legendre transformation F^i — F^i — (Jni^ni- Here, we use the sum convention for repeated indices 
for the sake of brevity. 

From now on, we restrict our considerations to linear isotropic elasticity. Taking the relaxed state of the a phase 
as reference, the elastic contribution to the free energy there reads 
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FIG. 1: Schematic figure of the steady state growth, with constant velocity u, of a thermodynamically favored phase /3 into 
a metastable phase a, which is initially at temperature Too. The motion of the interface is locally described by the interface- 
normal velocity which is given by the local heat conservation Eq. ([2]). The dashed lines at the two side borders indicate 
that we will discuss two different types of outer boundary conditions: Free growth, where the strains decay for x, y — > oo, and 
channel growth, with fixed displacements and thermal insulation at the walls. In the growth direction y the system is assumed 
to be infinite in both cases. 

where E and v are Young's elastic modulus and Poisson ratio, respectively. In contrast, the elastic contribution to 
the free energy density of the new phase /3 reads 
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which has an other state of zero elastic energy due to the lattice strain e^^, assigned to the phase transformation (see 
for example (author?) 0,0]). In the following we will call e"^. the eigenstrain of the transformation. Particular cases 
of eigenstrains and their interpretation will be discussed in the next section. 

The mechanical equilibrium and coherency conditions provide expressions for the discontinuous jumps of the strains 
at the interface 
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which are only related to the imposed eigenstrain. Defining an eigenstress tensor via Hooke's law 
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we obtain after a few straightforward algebraic manipulations for the elastic contribution to the local equilibrium 
condition Eq. ([3]) 
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Apparently, this expression depends on the strain state of the a phase at the interface. 

As already implicitly assumed above, the elastic degrees of freedom relax fast on diffusive timescales, and therefore 
the application of static elasticity is legitimate. Therefore, Newton's second law becomes 



dxk 



= 0. 
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FIG. 2: Visualisation of the hexagonal to orthorhombic transformation. From the symmetry its obvious, that the orientation 
angle § of the new orthorhombic phase can have only the three different values 0, ±27r/3. 

The above set of equations now has to be supplemented by boundary conditions at the external system boundaries, 
and we will discuss different scenarios below. They complete the self-sonsistent description of the moving boundary 
problem, which require the simultaneous solution of the elastic and diffusion equations, coupled via the boundary 
conditions at the propagating interfaces, leading to a complicated nonlinear and nonlocal problem. 

We note that we will later also discuss more complicated situations with more than one growing phase. Then, also 
similar boundary conditions and equations have to be set up for the interfaces between them. Typically, then also 
triple junctions will appear, where the contact angles are given by Young's law. 

III. MODELS OF STRUCTURAL TRANSITIONS 

The simplest case of a characteristic lattice strain is to assume the bond length of the new phase (3 to be uniformly 
longer or shorter in all directions in comparison to the reference phase, i.e. e'^j. = ef^. = eSik, which implies a density 
change. Notice, that this kind of transformations, to which we will refer to as dilatational eigenstrain, only changes 
the volume of an elementary cell, and not its relative dimensions or its scaled shape. 

The opposite case are transformations, where the volume of the unit cell is conserved and only the shape is changed. 
These shear transformations are characterized by a traceless eigenstrain tensor. A particular transition involving shear 
strain can occur e.g. in hexagonal crystals, where the transition leads to a lowering of the symmetry from C'e to C2. 
This is the case for example in hexagonal-orthorombic transitions in ferroelastics (see [l^ and references therein). 
Let the principal axis Cg be oriented in z direction. For simplicity we neglect all other possible strains with higher 
(axial) symmetry. By proper choice of the crystal orientation around the main axis in the initial phase, we obtain the 
new phase in three possible states in the xy plane due to the original hexagonal symmetry (See Fig. [2]) . Then the 
nonvanishing components of the eigenstrain tensor e°j. = eff. are = — ej^. = ecos2i? and e^y — esin2??, where the 
orientation angle d is either ?9 = or = ±27r/3. 

We also consider a more generic mixed mode case, where the self-strain of the transformation is a linear combination 
of the above two different transformations. We will discuss the following two mixing cases 

4 = 4('7) = ^4±(i-'y)4 (11) 

and refer to them as positive or negative mixing. 

It is important to mention that switching from positive to negative mixing corresponds to a rotation of the eigenstrain 
tensor by an angle tt/2. Since our model is isotropic in all respects except the choice of the eigenstrain tensor, the 
only predefined direction is given by the growth velocity. Therefore, we can interpret the negative mixing case to 
represent physically the same system as the positive mixing case but with a growth direction perpendicular to it. 
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IV. FREE GROWTH 



As a first application, we discuss the growth of the stable f3 phase from the metastable a phase in an infinitely large 
system, which is exposed to the temperature Tq^ far away from the interface. For simplicity, we discuss an effectively 
two-dimensional infinite system, assuming translational invariance in z direction. In particular, we assume a plane 
strain situation from point of view of elasticity. Also, we choose the elastic boundary conditions such that all stresses 
decay far away from the interface. 

Seeking for steady state solutions in free space, we use Green's function methods for solving both the elastic and 
the diffusional problem, to finally obtain self-consistently the shape and the velocityof the growing front. While the 
derivation of the integral representation of the thermal field is well-known, see e.g. IJ], we present the calculation of 
SFei in more detail [lol |. To this end, we express the elastic problem using the eigenstrain-independent stress tensor 

aik] for the a phase it is the same as the nomimal stress, crj^' = crjj!', whereas in the /? phase it is cr|^' = a\^} + a^f.. 
In terms of the new stress tensor dik, the mechanical equilibrium equation duik/dxk = ft introduces a force density 
which is localized on the interface and vanishes in the bulk, 

h = ^!f'--!:' = ^^- (12) 

This means that the problem of two coherently connected solids with different reference state is equivalent to a mono- 
lithic material with a distribution of point forces along the interface. We can therefore use an integral representation 
of the displacement field. 



J ds'G,k{r-r')fk{r'), (13) 



where the integration is performed along the interface, and the free-space Green's tensor Gik{r — r') of the elastic 
problem is given by (isj 

This function is the displacement response at r to a point force located in the origin. As a result, the mismatch 
provokes long-ranged elastic deformations, and therefore induces a nonlocal character to the problem. The above 
expressions can then be used to calculate the elastic influence on the interface temperature Wmt- 

Far behind the tip the interface profile becomes parabolic and is described by the asymptotic Ivantsov parabola 
[l^ . In this region, the elastic strains have decayed, and only the constant local contribution to SF^i (second term in 
Eq. [9]) remains. Consequently, the effective driving force A consists not only of the thermal undercooling A, which is 
controlled by the applied temperature Too, but also a contribution A^i from the elasticity, 

- 2(l-i.2)L2 • (15) 

We can consider A^i as a parameter quantifying the strength of the elastic effects. The elastic effects therefore induce 
a shift of the equilibrium temperature, and the effective undercooling A is lower than the actually applied thermal 
driving force, 

A = A-Aei, (16) 

which reflects the elastic hysteresis [ll B. 

By elimination of the thermal field we can derive the steady state equation for the shape of the solid-solid 
interface in a closed dimensionless representation, which reads in the co-moving frame of reference , 



dpK TeqCSFel _ p 
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dx' exp[-p{y{x) - y{x' j)]Ko{pfJ,{x, x')). (17) 



Here, p ~ vR/2D denotes the Peclet number {v is the steady state velocity), and it is related to the driving force 
by the implicit relation A = y^Trp erfc(p) . Furthermore, Kq is the modified Bessel function of third kind in zeroth 
order, and ^{x,x') = [{x — x'Y + {ui^) — 2/(2;'))^]^^^- the above expression we rescaled all lengths by the radius of 
curvature R of the asymptotic Ivantsov parabola. 

Consequently, the central equation ([T7]), together with the expression for the elastic contribution ©, where the 
strains are expressed via the integral representation (|13p results in a complicated integro-differential equation for 




FIG. 3: Phase field shapes of single /3 phases growing with constant velocity v into the metastable a phase in a finite channel 
of width W. The shapes in a) and b) correspond to the dilatational case {rj — 1), with A^i — 0.4 and an undercooling, which 
is chosen such that according to Eq. (|23|l it would correspond to an asymptotic phase fraction A = 0.65. While a) shows the 
shape of the symmetrical finger, 6) shows the shape of the sightly faster parity broken finger (see Fig. [5^). In panel c) we have 
plotted the shape in the case of negative mixing eigenstrain, ri — 0.5, which in that case corresponds to growth in the privileged 
direction. The parameters are Aei = 0.05 and A = 0.6 (for the resulting velocity see Fig. [SJa). 



the interface shape. It differs from conventional dendritic growth by the appearance of the elastic term. To better 
understand its influence, we briefly review the solvability conditions for classical dendritic growth [3, 0], for situations 
with pure dilatational and shear strain and cases with a triple junction In classical dendritic growth, any 
solution for isotropic surface tension with nonzero solvability parameter a = do /pR must exhibit a finite negative tip 
kink angle (f>. For this reason, dendritic growth is not possible for isotropic surface energy, since no smooth dendrite 
profiles exist. A dilatational eigenstrain simply shifts the equilibrium interface temperature, and therefore suffers 
from the same missing selection as the dendrite problem without surface tension anisotropy. The presence of a shear 
eigenstrain, which will also be discussed in more detail below, leads to the propagation of a bicrystal. In this case, 
the elastic effects operate as new selection mechanism. Independent of this elastically induced selection mechanisms, 
the presence of a triple junction alone also works as a selection mechanism. 

We solve the eigenvalue problem Eq. (|17p numerically for the steady state shape and the stability parameter a = a* 
in the spirit of Ref. [l^l , and the results are presented in Section IVII 



V. CHANNEL GROWTH 



The growth of the (3 phase in a channel of finite width W differs significantly from the free growth scenario discussed 
above. The specific behavior depends of course on the boundary conditions, and we assume fixed grip conditions on 
the long sides of the channel, i.e. the displacement is constant. Again, we discuss an effectively two-dimensional 
system by the assumption of translational invariance in z direction (see Fig. [3] for the setup). Furthermore, we assume 
that the a phase far in front of the growing tip is fully relaxed, which is in agreement with the boundary conditions 
Ux = and Uy = along the boundaries of the whole, infinitely long strip. From a thermal point of view, we assume 
thermal insulation along these boundaries, so no temperature flux through the surfaces, dT/dx = on the outer 
boundaries. Far ahead of the tip, the temperature is given by the dimensionless undercooling A. 

Due to the phase transition, the release of the latent heat leads to a heating of the system, therefore the temperature 
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far behind the tip T_oo is higher than in front. Under the given circumstances it is quite obvious that in this tail 
region a homogeneous state is reached. In particular, since the interfaces do not move there, the temperature becomes 
indeed constant. 

For simplicity, we discuss here only the specific case of a purely diagonal eigenstrain, i.e. eij = for i ^ j, to 
demonstrate the analytical procedure. In contrast, in the case of an eigenstrain tensor with nonzero off-diagonal 
elements, bicrystal patterns emerge, and we will briefly discuss this situation below; it is quite obvious that the 
analysis can be performed in an analogous way. 

Due to the confinement of the system, phase coexistence is possible in a finite range of undercoolings, and we can 
calculate the asymptotic volume fraction A of the new phase (see Fig. |3]). Notice, that in contrast to the isothermal 
situation discussed in 0], the temperature in the tail region is not a control parameter but has to be found self- 
consistently, due to the thermally insulating boundaries. In terms of the asymptotic phase fraction A, the width of 
the /? phase in the tail region is written as AM^. Mechanical compatibility requires 

(1 - A)eH + AeE! = 0, (18) 

since the average lateral strain must vanish due to the boundary conditions for the displacement. Together with 
the mechanical equilibrium condition of force balance, we can then calculate the elastic fields in the tail region, and 
consequently the elastic energy. Then, we calculate the energy excess bT ^ due to the appearance of a finite fraction 
of the new phase A. It is defined as the difference between the energy for A = without any phase /3 and the energy 
with an arbitrary but finite fraction, 

A 

6T^-W (\Ff + (1 - A)F[fl) + (^ ^(^e,-T(AO) ^ ^^^^ 

Initially, there is only the elastically relaxed but metastable phase a at the constant temperature A. As soon as we 
have a finite phase fraction A, the energy increases by elastic and capillary contributions. Finally, the integral appears 
due to thermal insulation, since an increase of the amount of the phase (3 by phase transformation processes generates 
latent heat, and therefore causes an increase of the temperature. Using the heat conservation condition we find the 
relation between this temperature and the phase fraction to be 

T{X)=T^ + X^. (20) 
Therefore, we obtain for the energy excess as a function of the asymptotic phase fraction A 

SH>^) = ^ ((A - AeOA - i (1 + ae^AeO - ^) , (21) 
where the parameter A^i, as already introduced above in Eq. (fTS]) . defines the strength of elastic effects, and 

_ CT,,(l + ^)(l-2^) KJ' 

L^{l-iy)E Aez ^ ' 

is a parameter depending on the type of eigenstrain. The maximum of this energy excess, dAF{X)/d\ — 0, provides 
the equilibrium fraction, which expresses the thermodynamic balance between elastic deformation and free energy 
release due to the phase transition, 



1 

1 + aelAel 



Aeg = , , „ ^ (A-AeO- (23) 



If the fraction Agg is in the range < Xeq < 1, we obtain phase coexistence in the asymptotic regime far behind 
the tip. We would like to point out that one can derive this equation also by using the phase equilibrium condition 
Eq. ©. 

Then growth of the f3 phase demands that the energy excess (f2T|) for the equilibrium fraction (|23p has to be positive, 

WL'^ /^l /, , _ A N x2 2do 



CTeq \ 2' ' W 



SH>^) = 7^ o (1 + «eiA,0 A^ - ^ > 0. (24) 



This is equivalent to the condition A > Xcrit, where Xcrit is given by 



,2 4do 1 



w l + ae,A,r (^^^ 
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We note that for the shear transformations discussed above, the most favorable process is the growth of a bicrystal, 
where two twinned phases (3 and /3' with d = ±2-1x1 "i grow together, because this is energetically most favorable Q- 
For the bicrystal configuration, the analysis can be performed in a similar way, and then also the (3/(3' grain boundary 
energy appears in the energy balance analogous to Eq. (|19p . 

Considering the channel geometry, we use a phase field formulation to solve the full free boundary problem numer- 
ically. This method is very versatile and allows to study even more complicated scenarios, e.g. the full dynamical 
evolution and not only the steady state regime, three-dimensionsal situations with full anisotropy of elasticity and 
surface tension or simulations beyond the symmetrical model, where the two phases have different diffusion constants. 
Nevertheless, for brevity of presentation, we confine our investigations to the simple, isotropic, two-dimensional and 
thermally symmetric case that we also used for the Green's function approach in the free space geometry and the 
above analytical calculations. 

The phase field method introduces an additional numerical lengthscale ^, the width of the transition region between 
the different solid phases. This means that an order parameter, which discriminates between the different phases, has 
a profile = [1 + tanh(n/^)]/2 in normal direction through an interface. We operate here in the sharp interface limit, 
which means that this lengthscale is significantly smaller than all other physical lengthscales present in the problem. 
Then the description recovers the desired sharp interface equations that were proposed above. 

Hence, we introduce the phase field having a value = 1 in the metastable initial phase, and a value cf) — Q in 
the new phase. We start from a free energy functional 

F[cj>, u^] ^ J{fs + fdu, + f) dV, (26) 

V 

where fs{'V4>) — 37^(V(/))^/2 is the gradient energy density and fdw{4') = 67(^)^(1 — 4>Y /£, is the double well potential, 
guaranteeing that the free energy functional has two local minima at = and 0=1 corresponding to the two 
distinct phases of the system. The form of the double well potential and the gradient energy density are chosen such 
that the phase field parameter ^ defines the interface width and the parameter 7 corresponds to the interface energy 
of the sharp interface description [l^. The free energy density f{4>,w,Ui) of our phase field description is chosen to 
be an interpolation between the two distinct free energies of the model, 

f{^,u,, w) ^F]^h{cp) + Ff{l - hm + L^— ^(1 - hm, (27) 

eq 

where we choose the interpolation function h to be h{4>) = 0^(3 — 20). It is the simplest polynom satisfying the 
necessary interpolation conditions h{Q) = and h[l) = 1, and having a vanishing slope at = and = 1, in order 
not to shift the bulk states. 

The evolution equation of the phase field is given by the variational expression 

90 M fSF\ 



37a ^''^ 

For large values of the mobility M we recover the desired case of diffusion limited growth. The elastodynamic equations 
read 

fSF\ da,k 

which recover static elasticity for slowly propagating interfaces in comparison to the sound speed Vs ^ {E/ p)^/"^ , with 
p being the mass density. 

For the temperature field we have the usual thermal diffusion equation with the motion of the phase field or interface 
being a source of latent heat, 

where the prime denotes the derivative with respect to 0. The phase field model presented here is very similar to the 
model in [l9|. However, there only the very special case of melting or crystallization of a confined sphere with elastic 
effects was discussed. 

For the simulations we typically use parameters W/do = 100 and W/(, — 80. 
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FIG. 4: The left panel shows the stability parameter a* (see Eq. (|3Hl ) as function of the mixing parameter -q for the privileged 
growth direction with an eigenstrain tensor when Aez/p = 1. On the right panel we show the stability parameter a*(Ae;/p) 
for two different values of mixing rj. 



VI. RESULTS 



Apart from the theoretical investigations above, we also performed numerical simulations for both the free growth 
in an infinite system and in a channel or narrow strip. As already mentioned before, we observe the growth of a single 
new phase if the eigenstrain tensor is purely diagonal, and twin structures for shear transformations with § = ±27r/3. 
Also, we consider mixed mode situations which are between these two extreme cases. 



A. Single crystal 

We first consider the situation of free growth. Here, the simplest case is that of a dilatational eigenstrain where 
all diagonal elements have the same value, i.e. e°j, = e5ik] this case is obtained from Eq. pT|) by setting r] — 1. As 
already mentioned before, selection of a steady solution in the absence of elastic effects is not possible for isotropic 
surface tension. It is quite remarkable that also the inclusion of elastic effects due to the dilatational eigenstrain of 
the transition does not change this situation; the reason is that the elastic effects can be mapped to a change of the 
undercooling only, and therefore we recover again the missing selection of conventional free dendritic growth ^10.] . It 
turns out that also for the volume-preserving shear strain transitions with e'^y = — e^^. = e, where all other components 
e^f. vanish, i.e. 77 = and = 0, steady state solution do not exist either. 

Under these circumstances it is very suprising, that for values of the mixing parameter in between zero and one 
0<7y<l,t? = 0wedo find steady state solutions for free growth and negative mixing. This is shown in Fig. |4^, 
where the selected "stability parameter" a = do/pR — a*, which is a dimensionless measure for the growth velocity, 
is plotted as function of the mixing parameter rj. In general, the eigenvalue a* depends on rj, A^i and p, but in the 
regime of small Peclet numbers the elastic driving force and p appear only in the combination Aei/p. Fig. [Dd shows 
the dependence of the selected stability parameter on A^i/p for two different values of 77. For the other type of mixing 
^tk corresponding to the unprivileged direction of growth, we do not find steady state solutions in free space. We note 
that from the selected stability parameter a* one can extract the steady state growth velocity v via the relation 

v = ^aV. (31) 

In contrast, for the growth in a channel, steady state solutions do also exist for the case of a dilatational eigenstrain, 
77 = 1 (see Fig. [SK)- As expressed through Eqs. ([24|) and (l25|) . the phase /3 can not grow for phase fractions A below 
the threshold Xcrit- Beyond this value the growth of a symmetrical finger (dendrite) is the most favorable solution. 
We note that in situations where free growth solutions do not exist, the still existing solutions in the channel must 
have a scale that is determined by the size W of the channel. For higher driving forces an asymmetric finger grows 
with higher velocities, and is therefore more likely to be observed; typical finger shapes are shown in Fig. [3l 

For the situation of 77 = 0.5, where steady state solutions exist also in the free space, we performed a quantitative 
comparison of the results of the two methods in Fig. For this comparison we directly relate the fraction A to the 
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FIG. 5: Channel growth with diagonal eigenstrain. The left panel shows the dependence of the growth velocity on the driving 
force A, where we choose a dilatational eigenstrain {rj = 1, A^i = 0.4). For small driving forces we obtain velocities and shapes 
corresponding to a symmetric solution (see Fig. |3ji) . For higher driving forces a first order transition towards asymmetrical 
solutions takes place (see parity broken finger morphology in Fig. [St). In the right panel the comparison of growth velocities 
between phase field and the free space calculations is shown, for the case of negative mixing corresponding to the privileged 
growth direction (ri = 0.5, i9 = 0, A^; — 0.05). The agreement is reasonably good, especially for driving forces A slightly above 
the critical phase fraction Xcrit of zero growth velocity. 

driving force A for free growth via Eq. ([23| and Eq. (fT6| . Then, in the limit of small Peclet numbers, the latter is 
related to A according to p w A^/tt. Finally, via the calculated function a*{Aei/p) and Eq. (|3ip we obtain the growth 
velocity as a function of A, which is then compared to the phase field results. 

The results for the growth velocity agree well in a range of intermediate fractions of the new phase /3. In particular, 
we confirmed that the steady state solutions which were obtained by the Green's function method are indeed stable, 
as the dynamical phase field method tracks only these stable branches. For low phase fractions, growth is not possible 
below the threshold Xcrit, which tends to zero for infinitely wide channels. For very high driving forces, the growing 
finger is influenced by the geometrical confinement of the channel width and also affected by kinetic corrections to 
the sharp interface equations of motion, which therefore lead to a deviation from the free growth result which was 
obtained by the Green's function method. 

It is interesting to compare the situation to a case without elasticity, where the problem reduces to classical growth 
in a channel [l^, [2l|, i23, 2A, 25] ■ It is known that there even in wide channels symmetric solutions do not exist below 
a fraction A = 1/2, whereas we found that the eigenstrain leads to the existence of solutions at already significantly 
lower driving forces A, see Fig. [5^. 

B. Bicrystal growth 

In contrast to the cases studied above, the presence of off-diagonal eigenstrain entries lead to the preferential growth 
of bicrystal patterns. This behavior can easily be understood intuitively, since the eigenstrain in the /3 phase favors 
the presence of off-diagonal strain. Then the boundary conditions of fixed displacements require that the surrounding 
a phase is sheared in the opposite direction. This implies a finite elastic energy which is stored far behind the tip. If, 
on the other hand, a second finger of the same size with an off-diagonal eigenstrain with opposite sign appears, the 
elastic energy contribution that is associated with the off-diagonal elements of the strain tensor vanishes. Consequently, 
the total elastic energy in the tail region is reduced, and the growth of the bicrystal therefore thermodynamically 
favored. 

An additional energy contribution emerges from this bicrystal configuration, since the boundary between (3 and /?' 
generates a grain boundary energy jgi, . We consider in particular the case of vanishingly small twin boundary energy, 
Igb ^ 7- This affects also the geometrical situation at the tip, which is now a triple junction, and the contact angle 
(j) becomes zero, so the overall finger has a smooth contour. Typical shapes of such bicrystal pattern are shown in 
Fig. [51 which are here related to the growth in a channel, but they look similar for free growth. We note that we 
simulated only half of the channel using the symmetry of the growing structures 0] . 

Considering the free space geometry, the first important observation is that in contrast to the single crystal scenario, 
steady state solution do exist also for pure shear, 77 = 0, which is already discussed in [10|]. New is here that we found 
the same morphology also in the channel geometry, showing surprisingly good agreement in the comparison of growth 
velocities, done in the same way as was done for the single crystal described above (see Fig. [7|). Again, we would like 
to point out that the stability of the bicrystal solution obtained from the Green's function technique is confirmed by 
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FIG. 6: Bicrystal patterns growing in a finite channel, a) shows the single tip solution and h) the twinned finger morphology, 
which corresponds to the most favorable solution branch for higher driving forces. The positive mixing chosen for a) and b) 
leads to growth along the privileged direction and the parameters are 77 = 0.1, Aei = 0.1, A = 0.6. The panel c) shows the shape 
of the slow growth direction for r; = 0.5, Agj = 0.2, A = 0.6. 



0.5 



0.4 



Q 0.3 
CvJ 

o 

p 0.2 



0.1 



— Free growth 
O Channel growth 
X Critical fraction A. 




FIG. 7: Bicrystal with pure shear eigenstrain, rj = G: Comparison of growth velocities obtained from the boundary integral 
technique for the free space geometry [10] with the velocities from phase field simulations for the channel geometry (Aej = 0.3). 



the phase field simulations. 

We also did investigations of the bicrystal free growth considering mixed eigenstrains, and the results are summarized 
in Fig. [8] and [9l In the first graph, the stability parameter tr*, which is again a measure for the dimensionless velocity, 
is shown as function of the admixture parameter 77 for fixed ratio Agj /p. As before, the elastic driving force and the 
p appear only the combination l^ei/p for small Peclet numbers. 

As already mentioned above it is quite interesting that the appearance of an additional dilatational strain leads to an 
effective anisotropy of the system: Whereas in the case of a positive admixture e+ the velocity grows with increasing 
•q, the opposite happens for negative admixture e~. This suggests that in free growth the bicrystal structure will 
eventually grow in the direction with the higher propagation velocity. For the fast branch, we did not find solutions 
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FIG. 8: Plot of the stability parameter o'*{r]) as function of the admixture parameter r], for the two different types of mixing 
and Aei/p = 1.0. The solid line corresponds to the privileged direction of growth, i.e. positive mixing, while the dotted line 
shows the dependence for the unprivileged negative mixing configuration. 




Ae|/p Ae|/p 

FIG. 9: The left panel shows the stability parameter a* as function of the control parameter A^i/p, for positive mixing, whereas 
the in the right panel a*{Aei/p) is plotted for negative mixing. Additionally, both graphs contain the (T*(Ae;/p) dependence 
for the pure shear configuration = 0, as a reference to compare both types of mixings, which was taken from [lOj]. 

beyond the point rj « 0.25. 

In Fig. [in] we show the comparison of the two geometries for the bicrystal case. Again, the phase field results 
confirm the stability of the solution branch for low driving forces. In this region, we obtain a convincing quantitative 
agreement between the growth velocities obtained for the free growth scenario using the Green's function method, 
and the channel geometry with the phase field method. In the left panel the favored direction is shown (positive 
admixture) for 77 = 0.1. In this case we again observe dendritic solutions for small driving forces A (see also Fig. 
For the chosen parameters in the channel geometry, all solutions on the upper branch are dendritic. On the lower 
branch, however, the steady state shapes are double- fingers, which means here the curvature directly at the point of 
symmetry has changed its sign and we obtain two tips (see Fig. Et"). The velocity of these branched solutions is here 
always smaller than for the single-tip structures. 

For the other direction, i.e. negative admixture, and a mixing parameter rj = 0.5, we show the comparison of growth 
velocities in the right panel of Fig. [TOl Here, we obtain a convincing agreement between the different methods and 
associated geometries at least for small driving forces, indicating dynamical stability also for this type of solution. 
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FIG. 10: Plot of the dimensionless growth velocity as function of the "driving force" A for bicrystal growth with mixed 
eigenstrain: comparison of free growth results obtained by Green's function methods to phase field simulations for the growth 
in a channel. Left panel: positive admixture with rj = 0.1. Right panel: negative admixture, rj = 0.5. 

Obviously, the growth velocities are here substantially lower than for the previous case, and for fully free growth an 
initial bicrystal seed would therefore preferentially grow in the "fast" direction. 

VII. SUMMARY AND CONCLUSION 

In summary, we investigated the diffusion-limited kinetics of solid-solid transformations in the presence of lattice 
strain. We studied in particular the free growth with Green's function methods and the growth in a channel by phase 
field techniques, exhibiting a comparable behavior. It turns out that the elastic effects have a strong influence on the 
selection of the arising patterns and the corresponding steady state velocities, which differ crucially from conventional 
dendritic growth. In particular, we find that the selected stability parameter and the corresponding growth velocity 
can be much higher compared to the classical case of dendritic growth, where the solutions are selected by the tiny 
effect of the anisotropy of surface tension. Different structural transitions, which are described by dilatational and 
shear eigenstrains, lead to a rich behavior, and we find very different structures already for the simple cases that were 
discussed here. 

Here, we focused on growth situations that are controlled by thermal diffusion, but it is known that the also 
interesting scenario of chemical diffusion can be mapped to this case. 

In comparison to the Green's function method, which is computationally very efficient, the phase field technique has 
the advantage of higher flexibility. In particular, it can easily be extended to three-dimensional cases, it can be used 
also for non-steady state situations or anisotropic or non-symmetrical situations with different material parameters 
in the phases. 

This work was supported by the DFG grant SPP 1296 and the German- Israeli Foundation. 
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